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ABSTRACT 

The detection of the high-frequency peaked BL Lac object (HBL) SHBL J001355.9-185406 (z=0.095) at high (HE; 100 MeV<E<300 GeV) and 
very high-energy (VHE; E > 100 GeV) with the Fermi Large Area Telescope (LAT) and the High Energy Stereoscopic System (H.E.S.S.) is 
reported. Dedicated observations were performed with the H.E.S.S. telescopes, leading to a detection at the 5.5 cr significance level. The measured 
flux above 310 GeV is (8.3 ± 1.7 stat ± 1.7 sys ) x 10~ 13 photons cirT 2 s -1 (about 0.6% of that of the Crab Nebula), and the power-law spectrum has 
a photon index of T = 3.4 ± 0.5 stat ± 0.2 sys . Using 3.5 years of publicly available Fermi-LXT data, a faint counterpart has been detected in the 
LAT data at the 5.5 cr significance level, with an integrated flux above 300 MeV of (9.3 ± 3.4 stat ± 0.8 sys ) x 10 -10 photons cm" 2 s and a photon 
index of T = 1.96 ± 0.20 stat ± 0.08 sys . X-ray observations with Swift-XRT allow the synchrotron peak energy in vF v representation to be located 
at ~ 1.0 keV. The broadband spectral energy distribution is modelled with a one-zone synchrotron self-Compton (SSC) model and the optical data 
by a black-body emission describing the thermal emission of the host galaxy. The derived parameters are typical of HBLs detected at VHE, with 
a particle-dominated jet. 
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1. Introduction 

Observations of blazars in y-rays offer the unique possibility 
to probe one of the most violent phenomena in the Universe. 
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Blazars are active galacti c nuclei (AGN) with th eir jets point- 
ing towards the observer (lUrrv & Padovanil 119951). Two classes 
have b een distinguished (see e.g. iGiommi et al Ghisellini 
1201 ll) : the fiat spectrum radio quasars (FSRQ) are the most pow- 
erful blazars and their optical spectra exhibit absorption and 
emission lines, whereas the BL Lacertae (BL Lac) class is less 
luminous and presents weaker emission lines, with the equiva- 
lent widths of the strong est ones smaller than 5A (IStickel et alj 
U99lHStockeetalJll991h . 

Since the discovery of the first extragalactic source in 
the VHE domain (Mrk 421, iPunch et al.fl 19921) . the num- 
ber of detected objects increased from a few in early 2000 
up to 50 objects at the beginning of 2013 (see TeVCafl for 
an up-to-date overview.). The advent of the current genera- 
tion of atmospheric Cherenkov telescopes (H.E.S.S., VERITAS, 
MAGIC) and the subsequent gain in sensitivity permit the de- 
tection of fainter and more distant sources. The efforts spent 
on unders tanding the properties of blazar s and predicting VHE 
emission dCostamante & Ghisellini 120021: iMassaro etalTl2008t 
lAbdo et aL 2010al) . together with measurements in the high- 
energy (HE; 100 MeV<E<300 GeV) range, greatly contribute 
to the detection of new sources. 

The spectral energy distribution (SED) of blazars is bimodal 
with one bump at low energy (from radio up to X-rays) and one 
at higher energy (from X-rays to TeV). BL Lac objects are fur- 
ther divided into two subclasses depending on the ratio betw een 
the X-ray and the radio fluxes dPadovani & Giommilll995l) : the 
low-frequency peaked BL Lac objects with a peak below UV 
wavelengths and the high-frequency peaked BL Lac (HBL) ob- 
jects for which the peak is in the UV or X-ray range. The latter 
subclass represents the bulk of the currently known extragalac- 
tic VHE y-ray emitters detected by atmospheric Cherenkov tele- 
scopes (34 out of 50 early 2013), an d almost 50% of the secon d 
Fermi catalogue of HE AGN (2LAC. lAckermann et al.ll201 id) . 

The low-energy bump of the SED is attributed to synchrotron 
emission of relativistic leptons (e + e~ pairs) moving along the 
jet. The origin of the high-energy component is less defi- 
nite. Leptonic models invoke inverse Compton (IC) scattering 
on ei ther the synchrotron ph otons (synchrotron self-Compton, 
e.g. lBand&Grindlavl[T98l or an external photo n field (ex- 
ternal Compton, e.g. Dermer & Schlickeiserlll993l) . The y-rays 
can also be produced by hadronic interac tions, such as photo- 
production of pions ( e.g. iMannheiml 19931) or synchrotron emis- 
sion of protons (e.g. I Aharonianll2000h . 

First detected in X-ra ys with ROSAT (1RXS J001356.6- 
185408. 1 Voges et alii 19961) . SHB L J001355.9-185406 w as later 
identified as a BL Lac object bv lSchwope et alj (20001) . With a 



radio flux of 29.6 mJy at 1.4 GHz (ICondon et al 



19981) and an 



X-ray flux, betwee n 0.1-2.4 keV, of 1.26 x 10~"erg cnT 2 s -1 
(I Voges et all 1996b . it fulfilled the selection criteria of the seden- 
tary survey of extre me high-energy peaked BL Lacs (SHBL, 
IGiommi et al.ll200"5l) . making this source a member of the HBL 

subclass. 

Its relative proximity (z=0.095. lJones et"al 12009b and its X- 
ray and radio fluxes are criteria that make SHBL J001355.9- 
185406 an interesting target for very high-energy (VHE , 
E > 100 GeV) observations dCostamante & Ghisellrnl |2002|) . 
Consequently, the source was observed with the High Energy 
Stereoscopic System (H.E.S.S.) and indeed repo rted as a TeV 
y-ray emitter in November 2010 dHofma nn 2010). 

The Fermi Large Area Telescope (LAT), launched 
on June 11, 2008, did not detect a high-energy (HE; 



100 MeV<E<300 GeV) counterpart after two years of op - 
eration (Fermi two-year catalogue, 2FGL, iNolan et all 120121) . 
However, the analysis of 3.5 years of data reported here reveals 
the presence of a faint source that is positionally coincident with 
SHBL J001355.9-185406. 

The H.E.S.S. and Fermi-LAT data analyses and results are 
presented in sections |2~T1 and l2~2l The multi-wavelength data set 
is presented in sections |231 and l2~4l for the X-ray and UV obser- 
vations with Swift and in section [231 for the optical observation 
with ATOM. The discussion in section[3]focusses on the descrip- 
tion of the SED in the framework of a synchrotron self-Compton 
(SSC) model. 

Throughout this paper a ACDM cosmology with Ho — 71 km 
s Mpc , Qa = 0.73 and Qm = 0.27 is assumed, resultin g in a 
luminosity distance of D L = 431 Mpc dHogg. D. W. Ill999t ?). 



1 http://tevcat.uchicago.edu 



2. Observations and analyses 

2.1. H.E.S.S. data set and analysis 

H.E.S.S. is located in the Khomas Highland, Namibia 
(23° 16' 18" S, 16°30'01" E), at an altitude of 1800 m above sea 
level. H.E.S.S. is an array of imaging atmospheric Cherenkov 
telescopes. Each of four telescopes used in this study (H.E.S.S. 
phase 1 ) consists of a segmented 13 m diameter optical reflector 
ilBernlohr et~aT1l2003h and a camera composed of 960 photomul- 
tipliers. The system works in a coincidence mode re quiring the 
detec tion of an air shower by at least two telescopes dFunk et al.l 
12004 . 

SHBL J001355.9-185406 was observed with H.E.S.S. be- 
tween MJD 54653 and MJD 55912 (July, 6 2008 - December, 
17 2011). Data were selected using the stand ard quality crite- 
ria (g ood weather, no hardware problem, see lAharonian et al.l 
120061) . yielding an exposure of 41.5 h acceptance corrected 
live time at a mean zenith ang le of 12.9°. The Model analy- 
sis dde Naurois & Rollandll2009l) was performed with the stan- 
dard cuts (threshold of 60 photo-electrons), leading to an en- 
ergy threshold of E± = 3 10 GeV for this observation. An excess 
of 153 y-ray candidat es has been found usi ng the Reflected- 
background method dAharonian et al.1 120061) to subtract the 
background. The total numbers of ON- and OFF-source events 
are Non = 830 and Noff = 8190, respectively, with a back- 
ground normalization factor a ^ 0.08 3. The source i s detected 
with a significance of 5.5 cr (Eq. 17 of iLi & Malll983l) . 

The distribution of excess events in a 2D map of 2° field of 
view centred on the source coordinates is given in Fig. Q] The fit 
of a point-like source, convolved with the H.E.S.S. point spread 
function (PSF), to the data results in the best-fit position of 
ffnooo = 00 h 13 m 52 s ±1.5^ + 1.3| ys and <5 J20 oo = -18° 53' 29" ± 
22" stat ± 20" sys , systematic uncertainties arising from telescope 
pointing. This is less than 2cr stat away from the test position 
of SHBL J001355.9-185406 , derived from radio observations 
( anm) = QQ h 13 m 56 s + 0.05 s and dnmo = -18° 54' 06" + 0.7", 
ICondon et al1fl998h . The source is compatible with a point-like 
source within the systematic uncertainties on the H.E.S.S. PSF. 
The distribution of ON-source events as a function of the square 
of the angular distance to the radio position is shown in Fig. [2] 
along with the OFF-source distribution. 
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The VHE spectrum (Fig. [3}, der ived from the H.E. S.S. data 
using the forward-folding method in lPiron et al 1 (12001b . is com- 
patibly with a simple power law of the form 



dN 
dE 



(1.16 + 0.23 st a,±0.23 sys )x 10 



-12 



-3.4±0.5 !tat ±0.2 s y ! 



cm 



where the decorrelation energy is 
grated flux above 310 GeV is (8.3 ± 



= 5 10 GeV. The inte- 

1.7 stat ± 1.7 sys ) x 10- 13 



photon s cm 2 s 1 . The 1 cr error contour was computed using 
Eq. 1 in lAbdo et al 1(l2010bh. and the systematic errors were eval- 
uated following lAharonian et al.l (|2006|). The analysis was c ross- 
checked with an independent method (IBecherini et al.l201lT) giv- 
ing compatible results. 

A search for variability was performed on a periodQ time 
scale. The corresponding light-curve (Fig. |4] top panel) does not 
show any significant variations exceeding the experimental un- 
certainties, and the fit with a constant function yie lds a y 2 proba- 
bility of 7%. The fractional excess variance F var dVaughan et al.l 
120031) is compatible with zero at the 2 cr level, and the 99% con- 
fidence level upper-limit is F val < 2.57. 
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Fig. 2. Distribution of ON-source events (circles) as a function 
of the square of the angular distance 9 2 from the source posi- 
tion. A cut at 8 = 0.1° is used to define the ON region with 
standard cuts. The distribution of normalized OFF-source (solid 
histogram) is flat and matches the distribution of the ON-source 
events at large 2 , which reflects a proper background subtrac- 
tion. 
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Fig. 1. Map of the y-ray excess measured with the H.E.S.S. tele- 
scopes around the position of SHBL J001355. 9-185406 in right 
ascension and declination (J2000). The map is smoothed with 
the H.E.S.S. PSF, shown in inset. 



2.2. Fermi-LAT data set and analysis 

The LAT on-board the Fermi satellite is a pair conversion tele- 
scope sensitive to y-rays between 20 MeV and 300 GeV. Its main 
charact eristics and performance can be found in lAtwood et al.l 
(I2009ah . The bulk of LAT observations are performed in an all- 
sky survey mode, where all objects are seen for about 30 minutes 

every 3 hours. 

Events passing the SOURCE selection (lAckermann et al.l 

120 12l) with a reconstructed energy between 300 MeV and 
300 GeV have been considered in this analysis. The correspond- 
ing instrumental response functions (IRFs) P7SOURCE_V6 and 



The^ 2 /NDF is 25.1/19 for a corresponding probability of 15%. 
A period is defined as the time between two full moons (lunation). 
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Fig. 3. Spectrum of SHBL J001355.9-185406 measured with 
H.E.S.S. using the Model analysis. Top panel: The butterfly rep- 
resents the 1 cr contour for the best-fit mod el. The spectrum i s 
obtained with the forward-folding method in lPiron et al 1 (120011) . 
The data points are derived a posteriori and should be consid- 
ered as residuals. For bins with a significance below 2cr, upper 
limits at the 95% confidence level are computed. The inset gives 
the 1, 2, and 3 cr confidence levels in the power-law index vs dif- 
ferential flux at 1 TeV plane. Bottom panel: Residuals of the fit, 
with the binning used in the forward-folding method. 



the latest public version of the ScienceTools (v9r27pl), avail- 
able from the Fermi Science Support CenteiQ web site, were 
used. A region of interest (ROI) of 15° radius around the ra- 
dio coordinates of SHBL J001355.9-185406 was defined to per- 



http : //f ermi . gsf c . nasa . gov/ssc/data/ analysis/ 
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form a binned analysis dMattox et al.lll996h . implemented in the 
gtlike tool. Additionally, cuts were applied on the rocking an- 
gle of the spacecraft, which was required to be smaller than 52°, 
and on the zenith angle of the events, required to be smaller than 
100°. 

The model used to describe the observed emission consists 
of all the sources contained in the 2FGL catalogue within this 
ROI, an isotropic extragalactic diffuse component described with 
the file iso_p7v6source . txt, and the standard Galactic model 
gal_2yearp7v6_vQ. fits. At the position of SHBL J001355.9- 
185406, a point-like source with a power-law model was added. 
The model of SHBL J001355.9-185406 has two free parame- 
ters: the differential flux at a fixed energy and the photon in- 
dex T. The spectral parameters of sources within 3° of the 
SHBL J001355.9-185406 position, as well as the diffuse back- 
ground normalizations and the variable and/or bright sources 
(namely 2FGL J2330.9-2144 and 2FGL J2345.0-1553), were 
left free to vary during the fitting procedure, whereas the other 
parameters were frozen to the 2FGL values. 

Using 1290 days (^ 3.5 years) of data, from August, 
4 2008 to February, 12 2012 (MJD 54682 to MJD 55973), 
SHBL J001355.9-185406 is detected with a test statistical (TS) 
of 34, which approximately corresponds to 5.5 cr. This source 
was not detected by Fermi-LAT at the time of the 2FGL cata- 
logue, which can be explained by its faint y-ray flux. The coor- 
dinates derived with gtf indsrc are compatible with the radio 
position within 1 cr. The source spectrum is described well by 
a power law with a differential flux of (3.5 + 0.9 stat + 0.7 sys ) x 
10~ 14 crrT 2 s~' MeV~' at the decorrelation energy E&x = 2.90 
GeV. The photon index of the source is found to be T = 
1.96 ± 0.20 stat ± 0.08 sys , yielding a total flux above 300 MeV 
of (9.3 ± 3.4 stat ± 0.8 sys ) x 10~ 10 photons cm" 2 s" 1 . The photon 
with the highest detected energy that is very likely associated 
with the source^ has an energy E » 37 GeV. Systematic uncer- 
tainties have been e valuated using the IRFs bracketing method 
(lAbdo et al.l 12009b). A fit with a log-parabola function does 
not improve the overall results. The differential flux has also 
been measured in four energy bins (see Fig. [3} by performing 
a gtlike analysis for which the photon index of the source is 
frozen to the best-fit value. The last bin has a TS<8 and a 95% 
confidence level upper limit has been computed. 

Due to the faint emission of the source in the 300 MeV- 
300 GeV range, a light-curve with six-months-wide time-bins 
was computed using the gtlike analysis chain. This light curve 
is presented in the second panel of Fig. [4] Upper limits at the 
95% confidence level were calculated for time bins with TS<4. 
No significant variability can be measured on this time scale with 
F var < 1.83. 

2.3. Swift-XRT data set and analysis 

X-ray observations are usually an important part of characteriz- 
ing the SED of blazars, since they probe the synchrotron radia- 
tion of the highest energy leptons. Target-of-opportunity obser- 
vations with the space-based Swift dBurrows et al.luOOl X-ray 
observatory have been conducted in September 2010. Data taken 
in photon-counting (PC) mode are processed with the standard 
xrtpipeline tool from the HEASOFT V6. 12, where a King func- 



tion fit to the PSF shows no evidence of any pile-up in any of the 
four observations. Events and background extraction regions are 
defined with a 60-pixel radius circle (corresponding to » 142"), 
with the latter centred near the former without overlapping . The 
Swift-XRT spectrum has been rebinned so as to have at least 
20 counts per bin using grppha, yielding a usable energy range 
between 0.3 and 9.0 keV for the summed spectrum ( 0.3 and 
5 - 7 keV for single observations). The weighted average column 
density of Galactic HI Nh = 2.13 x 10 20 cirT 2 has been extracted 
from the Leiden/Argentine/Bonn (LAB) Survey dKalberla et al.l 
12001 using the nH tool from HEASARCQ Multiple model spec- 
tra are tested with pyXspec 1 . 0, using the response functions 
swxpc(8tol2s6_20QlQlQlv013 and a dedicated ancillary re- 
sponse function (ARF) using xrtmkarf within FTOOLS at the 
location of the source in the field of view and with the summed 
exposure maps of the single observations. A single power law 
F(E) = KE r " poorly represents the fitted summed spectra , 
with^ 2 = 1.23(179). An F-test probability (lBevingtoiJ[l969l) 
of 5 x 10~ 4 prefers the simplest smooth curved function, a three- 
parameter log -parabola F(E) = KE~ a ~ blo&( - K> , for which all pa- 
rameters are given in Table Q] including the spectrum obtained 
after summing all four event files in xselect V2 . 4b and build- 
ing an exposure map with ximage V4 .5.1. Similar conclusions 
were drawn by iMassaro et al.l (1201 ll) on this object with a dif- 
ferent calibration. The source shows no indication of variability 
over the span of eight days (third panel of Fig. 

2.4. Swift-UVOT data set and analysis 

The UVOT instrument dRoming et al.l200l on-board Swift mea- 
sured the UV emission of SHBL J001355.9-185406in the bands 
U (345 nm) and UVW2 (188 nm) simultaneously with the X-ray 
telescope with an exposure of ~ 1.5 - 2 ks each. The instru- 
mental ma gnitudes and the co rresponding fluxes (for conversion 
factors, see lPoole et al.ll2008b are calculated with uvotmaghist 
taking all photons into account from a circular region with a ra- 
dius of 5" (standard aperture for all filters). An appropriate back- 
ground was determined from a circular region with radius 40" 
near the source region without contamination of neighbouring 
sources. The measured UV fluxes are corrected for dust absorp- 
tion using E(B - V) = 0.0246 magB (Schlege l et al. 1998) and 
the A A /E(B - V) ratios given in lSeatonl(ll979l) . 

The flux (Fig. [4] fourth panel), binned per observations, does 
not show any sign of variability in either of the two bands. 
Table [2] gives the chance probability obtained when fitting the 
data with a constant, as well as the average flux obtained when 
analysing all the observations together. 

2.5. ATOM data set and analysis 

The 75 cm Automa tic Telescope for Optical Monitoring (ATOM, 
IHauser et alj|2004l) . located on the H.E.S.S. site has been used 
to monitor the optical emission of SHBL J00 1355.9-185406 in 
Bessell B, R, and I filter bands (lBesselllll99dh over the last five 
years. The presented data have been obtained from MJD 54629 
to MJD 55897 (June, 6 2008 - December, 2 2011). A total of 
138, 188, and 2 observations in B, R, and I bands, respectively, 
were carried out with an aperture of 4" radius. The data were 



corrected for Galactic absorption using A^B) - 0.107 mag, 

5 see |http : //f ermi . gsfc . nasa . gov/ssc/data/analy sis/docume&oaffiloff/di<^^ WstfkkCEe^HWew . htm] 

for a definition. 1998) and the Bessell zero points (Bessell 1990) are used to 



6 A photon is considered as associated to the source if its recon- 
structed direction is within the 95% containment radius of the Fermi- 
LAT PSF. 



http : //heasarc . gsf c . nasa . gov/cgi-bin/Tools/w3nh/w3nh . pi 
|http : // irsa . ipac . caltech . edu / applications/DUST/| 
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Table 1. Parameters of power-law or log-parabolic fits to all four Swift-XRT observations, as well as the summed spectrum. Column 
4 gives either the power-law index Fx or the log-parabola value of a depending on which model best fits the data. Column 7 shows 
the reduced^ 2 value ;^ 2 , column 8 the F-test probability for the log-parabolic model when it is preferred over the power law, and 
column 9 the estimation for the unabsorbed flux (using cf lux) in the 0.3 - 10 keV range. The values of b and F-test are not provided 
when the log-parabola is not preferred. 



ObsID 


Start time 


exposure 


T x 


or a 


b 


K 




F-test 


^0.3-10keV 




MJD-55000 


M 








[IQ- 3 cm- 2 s-'keV- 1 ] 






[10-" ergcrrr 2 s~'] 


00031806002 


449.05 


3783 


1.95 


±0.06 


30 +U14 

V - J -0.13 


2.62 ± 0.09 


1.15(55) 


4 x 10 2 


1.28 ± 0.06 


00031806003 


449.98 


4190 


2.03 


±0.04 


2.46 ± 0.07 


1.28(61) 




1.38 ± 0.04 


00031806004 


451.05 


3848 


1.98 


±0.05 


26 +014 

—0. ! 3 


2.46 ± 0.09 


0.98(52) 


5 x 10~ 2 


1.30 ±0.08 


00031806005 


451.58 


2745 


2.02 


±0.05 


2.38 ± 0.08 


0.80(39) 




1.33 ± 0.09 


sum 




14566 


2.00 


±0.03 


0.23 ± 0.06 


2.73 ± 0.05 


1.15(178) 


5 x 10- 4 


1.36 ±0.03 



Table 2. Swift-UVOT U and UVW2 measurements and time- 
averaged optical magnitudes and fluxes measured with ATOM in 
R, B, and I Bessel filters. The last column gives the probability 
obtained when fitting the time series with a constant function. 
Magnitudes are not corrected for Galact ic extinction, wherea s 
fluxes are dereddened using the values o f lSchlegel et all d 19981) . 



Filter 


Energy 


magnitude 


Flux 






[eV] 




[10- 12 erg cm" 2 s" 1 ] 


% 


UVW2 


6.6 


17.4 


1.99 


63 


U 


3.6 


17.2 


1.44 


94 


B 


2.88 


17.5 


3.26 


99 


R 


1.77 


15.8 


6.09 


19 


I 


1.37 


14.2 


8.36 





convert the magnitude into flux units. The time-averaged flux 
in each band is given in Table [2] along with the corresponding 
energy. 

Table [2] also gives the % 2 probability obtained when fitting 
the time series with a constant function. The flux of the source 
(Fig- III fifth panel) is compatible with being constant over time 
in each band. 



3. Discussion 

The broadband data presented in this work from ATOM, Swift, 
Fermi, and H.E.S.S., together wi th archival data from th e NVSS 
Condon et all [19981) . 2MASS (ISkrutskie et all I2006T) . 6dFGs 



Jone s et al. " 120091) . and RASS (IVoges et a U 119961) catalogues 



are used to constrain physical characteristics of the source. 
Contemporaneou s data from the W ide-field Infrared Survey 
Explorer {WISE, IWright et al.1l2010l) in the bands 3.4, 4.6, and 
12 //m are also presented. The SED exhibits the usual two com- 
ponents, as can be noticed from Fig. [5] 

The Swift-XRT data allow the synchrotron peak position to 
be located at E StPe ^ - 1.00 + 0.07 keV, a quite common fea- 
ture of many o bjects in the SHBL cat alogue, that peak in X-ray 
as mentioned bv lGiommi et al.1d2005l) . Unfortunately, the source 
is too fain t for the ene r gy of t he IC peak to be measured accu- 
rately, but lAbdo et akl d2010al) found an empirical relation be- 
tween the Fermi photon index and the peak position of the IC 
component. Applying the relation to SHBL J001355.9-185406 
leads to £ ic , peak = 3.1^ Ge\0. 

The spectral index found between the optical and UV wave- 
band qt -uv is - _ 2.3 which is not in agreement with the one 



9 The use of a relation derived on a population of sources might lead 
to large uncertainties but in the case of SHBL J001355.9-185406 , a 
factor 10 in £i Cipe ak would not affect the conclusions of this work. 



observed between UV and X-ray, a uv - x - 1.5. The negative 
value of the former indicates that it does not originate from 
synchrotron emission of non-thermal particles but rather from 
the host galaxy thermal emission. This is supported by the non- 
variability of the ATOM light-curves and the good agreement of 
the optical flux with the archival data. Pl acing the source in a 
[3.4]-[4.6]-[12]yum colour-colour diagram dMassaro et al.ll20TT1) 
reveals that the infrared measurement is also dominated by this 
thermal emission. The host galaxy emissio n is modelled using a 
black-body, as in iKatarzvhski et al.l d2003l) . with a temperature 
of 4500 K and a total luminosity of 1.43 X 10 44 ergs 1 , which 
corresponds to 3.7 x 10 10 times the luminosity of the Sun. 

The non-thermal emission is modelled with a sim- 
ple, time-independent, one-zone homogeneous SSC model 
(lBand&Grindlavlll985l) . The emission zone is spherical with 
a radius R, moving relativistically towards the Earth with a 
Doppler factor 5. A population of leptons with a density Ni(y) 
is responsible for the synchrotron emission by interacting with a 
uniform magnetic field B and for the y-ray emission by IC scat- 
tering of the synchrotron photons. 

Some properties of the SSC model can be assessed using the 
available data. A constraint on the maximal energy of the leptons 
can be found assuming the last bin in energy of Swift-XRT is a 
lower limit for the maximu m energy reached by the sy nchrotron 
process E sm . dx . This gives dRvbicki & Lightmanll 19791) 



15(1 + z) V3j 



> 8keV, 



(1) 



m e 



where q is the charge of the electron and B is in Gauss. Equation 
□implies y max > 9.7 x 10 4 fi- 1/2 (T 1/2 . 

Since the same population of particles is responsible for 
the two spectral components and under the hypothesis that 
the IC scattering occ urs in the Thomson regime, Eq. 4 of 
iTavecchio et al.l d 19981) reads as 



8.6 x 10 7 £ 2 , 
B6 = (l+ z) — = 3.1 x 10 4 G, 

^ic,peak 



(2) 



where the energies are expressed in eV. For a magnetic field 
B of 0.1 G, this leads to a Doppler factor of 3.1 x 10 5 , which 
is unrealistic. For the photons of energy higher than Sic, peak, 
the scattering must occur in the Klein-Nishina regime in which 
case the leptons producing the highest energy detected follow 
£ic,max = ^Tmax'«eC 2 = 18TeV, which with the above con- 
straint on y m ax, in turn yields B6~ l > 6.3 x 10~ 8 G. 

Even with the previous calculations, the SSC model is still 
not fully constrained by the data. To reduce the number of pa- 
rameters, the density of leptons is described by a power law with 
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Table 3. SSC parameters used to reproduce the SED of 
SHBL J001355.9-185406. The calculation is shown as a black 
line in Fig. [5] 



Parameters 


Value 


B[G] 


0.05 


R [cm] 


3.5 x 10 16 


S 


10 


P 


2.2 


ymin 


1 


Tcul 


5.0 x 10 5 


N m 


2.0 x 10 53 



an exponential cut-off of the form N e {y) oc y~P ■ exp(-y/y cut ) 
which only has three parameters. An equipartition factor Q (de- 
fined as the ratio between the density of particle kinetic en- 
ergy and the magnetic energy density (i.e. Q = u,,/ub)), that 
is as close to unity as possible has also been require d, to mini- 
mize the total energy of the jet (see lGhisellinill2012h . With this 
model, the minimum factor allowing a fit to the data is Q = 50, 
which rules out the equipartition and the jet is found to be 
particle-dominated. The same c onclusions have bee n drawn for 
other sources, such as Mrk 421 (lAbdo et al.ll201 lbb or Mrk 501 
(lAbdo et alJliOlTah . 

The SSC calculation, together with the thermal black- 
body component, is shown in Fig. [5] and parameters are 
given in Table [3] The attenuation by the extragalactic back- 
ground light has been ta ken into account using the model of 
iFranceschini et al.1 (120081) . The index of the leptonic distribution 
is found to be p — 2.2 with a break at y cut = 5.0 x 10 , for a total 
number of leptons of 2.0 x 10 53 . The size of the emission region 
is R = 3.5 x 10 16 cm, the Doppler factor 5 = 10, and the magnetic 
field B = 0.05 G (i.e. Ub = 10~ 4 erg cm -3 ), values that comply 
with the above-mentioned constraints. The minimum variability 
time scale achievable within this scenario is 1.4 days, which can- 
not be tested with H.E.S.S. or Fermi given the flux of the source. 

To check that y cut can be related to the synchrotron cooling 
time or be acceleration effects, the comparison between the ra- 
tio of the light cross ing time t c in the e mission zone rest-frame 
and the cooling timq^j r coo i can be used. iTavecchio et al.l ([1998) 
suggest that this ratio should be between 1 and 3. For the model 
presented here, this ratio is close to 1.8 and so is not possible to 
disentangle a radiative cooling break from a cut-off in the lepton 
distribution due to acceleration effects. 



4. Conclusions 

Dedicated observations using the H.E.S.S. telescopes have re- 
vealed a new HBL, SHBL J001355.9-185406, as a VHE y-ray 
emitter, with a significance of 5.5 cr. The source has a flux of 
about 0.6% of the Crab Nebula flux above 310 GeV, with a soft 
photon index F = 3.4 + 0.5 sta t ± 0.2 sys . 

Using 3.5 years of Fermi data, the presence of a previously 
undetected counterpart in the HE range was found. Both HE and 
VHE spectra connect smoothly, leading to the conclusion that 
the same population of particles is likely to be responsible for 
the y-ray emission from 300 MeV to a few TeV. 

At lower energy, the optical measurements, which are con- 
temporaneous with the H.E.S.S. and Fermi-LAT observations, 
are found to be constant over time and compatible with the data 



10 



- LAT { 0.3-300 GeV ) 
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ATOM R band 



ATOM B band 



1000 1200 1400 

Time (MJD - 54000) 



The inverse of the cooling time for synchrotron dominated models 



is defined as t 



4 (TJC 

3 m„c- ' 



Fig. 4. Light-curves from observations at different wavelengths. 
From top to bottom: ;) H.E.S.S. period-by-period integral-flux 
above 310 GeV. The dashed line is the mean flux level. The ver- 
tical bars give the statistical errors, it) Fermi-LAT flux in the 
300 MeV-300 GeV energy range. For time bins where the test- 
statistic is below 4, an upper limit at 95% confidence level is 
reported. Hi) Swift-XRT flux in the 0.3-10 keV band, iv) Swift- 
UVOT measurement in U and UVW2 bands, v) ATOM magni- 
tudes in the B and R bands. The insets of panel Hi and iv present 
a zoom in time around the Swift observations. 

taken over the past decade. They are interpreted as thermal emis- 
sion of the host galaxy and successfully reproduced by a black- 
body model. 

An SSC model has been used to reproduce the non-thermal 
emission from radio to TeV energies. In this model, the y-rays 
with energies above the y-ray peak are produced by IC scattering 
showing evidence for Klein-Nishina suppression. Equipartition 
between the kinetic and the magnetic energy density is ruled out 
by the calculation, implying that the jet is dominated by the lep- 
ton kinetic energy. 

The number of extragalactic sources jointly detected by the 
current generation of atmospheric Cherenkov telescopes and the 
Fermi-LAT is increasing rapidly, allowing the y-ray sky from 
100 MeV to a few TeV to be probed and the mechanisms re- 
sponsible for the electromagnetic emission to be constrained. 
The knowledge of the non-thermal sky will increase with the 
advent of H.E.S.S. 2 and the future Cherenkov Telescope Array 
(CTA). The ten-fold increased sensitivity of CTA with respect 
to current generation atmospheric Cherenkov telesc opes and the 
possib ility of performing an extragalactic survey dDubus et al.l 
2012) will allow detection of hundreds of sources with fluxes of 
1 % of that the Crab Nebula. Such survey, together with ten years 
of Fermi-LAT data, will allow detailed population studies. 
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E [MeV] 

Fig. 5. Radio-to-TeV spectral energy distribution of 
SHBL J001355.9-185406 with ATOM (red circles, Table 0, 
Swift-WOT (green circles), Swift -XRT (green cross + brown 
lines, Table [T]), Fermi-LAT (blue), and H.E.S.S. (red) measure- 
ments. Archival data from NVSS, 2MASS, 6dFGS, and RASS 
are shown in grey. The dashed line is the black-body spectrum 
from the host galaxy at a temperature of 4500 K, and the black 
line shows the sum of the SSC calculation and the black-body 
spectrum. 
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